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Abstract 



Two photon decay widths of the J p = + scalar mesons a (980), / (980), / (1370) and 
Xco are calculated in a covariant model which is characterized by the quark - antiquark 
structure. Previously such models were used to calculate current form factors. Here a 



^ ■ different application is tried. A simple version of the model uses adjusted nonrelativistic 
model parameters with small quark masses. The results seem to prefer nonideal mixing 
of /q(980) and /o(1370). The calculated decay rate of Xco agrees with the experimental 
results. 
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1 Introduction 



The scalar mesons appearing around 1 GeV mass scale seems to be least understood 
among spin zero mesons. The experimental dataQ, on strong decays of a (980), / (980) 
and /o(1370) do not lead to the final understanding of the real structure of these mesons 
P| — Hl7|. Ideas exists that these states are KK molecules [@] — @. There is a suggestion 



that cio (980) has a four-quark structure with the strange quarks (ss) contribution [|16 



based on the analyses of the — > 7ao(980) — > 77777 flPT|. Such uncertainties suggest fur- 
ther investigation, which were concentrated on 27 decays of a (980), / (980), / (1370) 
and Xco mesons. The aim was to explain the experimental decay widths and at the 
same time to reproduce the measured masses. 

The relativistic quark model is used to correlate various data and to establish the 
connection between masses and decay widths. 

For that purpose one employs the covariant model |I| [19|] which includes the heavy 
- quark symmetry. As shown in the following section the model and the calculations are 



both covariant and gauge invariant. This model is a covariant generalization [[18], [19| 



of the well known ISGW model pOfl . However, here the usage of small quark masses 
is investigated, which means the avoidance of the weak binding limit approximation in 
its strictest sense ||20|| . That might better mimick the real quark fields which should 
appear in the photon emitting quark loop (Fig. 1) in the first order of QED/QCD 
expansion. In a very simplified version of that model, which is employed here, only 
the quark momentum distribution parameter j3 and the model quark masses appear. 
All parameters are correlated and compared with the nonrelativistic choices [EDI fTD, 122 



The quark masses, were treated as fitting parameters in a limited sense. Suitable values, 
allowed within the experimental uncertainty in current quark masses [[!]], were selected. 
That parameterization is expected to lead to a reasonable reproduction of meson masses. 
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No additional fitting was allowed when widths were calculated. In that way one can 
reproduce the measured T(27) reasonably well and make the prediction for the / (1370) 
27 decay. 

All results are based on the valent quark qq structure, which is characteristic of the 
model. The importance of qq structure has often been mentioned J7| — [rj|. Our basic 
loop diagrams, Fig. 1 below, correspond closely to the quark loop diagrams shown 
in Fig. 1 of Ref.(lO). Thus it is not surprising that our results, following from the 
somewhat more complicated diagrams, depends strongly on quark masses. 

Moreover our model contains the sea component constrained by the requirement of 



general Lorentz covariance, valid in an frame [18[ |19 |. In that way, indirectly, some 



other QCD structures, as discussed earlier [18|, [19[ , enter into our description. 

Predictions based on our simplified model version test how well, or how badly the 
model mimicks the real QED/QCD world. The quarkonium approximation [ f20f 



is 



investigated here in the circumstances which are different from the usual form factors 
related problems [p~8|, |2H|| . 

The results depend also on the quark flavor structure of the scalar mesons. They 



do distinguish among various propose uu, dd and ss mixing [D5j, |TJJ in f$ mesons. We 
investigate only lower laying states without entering into discussion of the states as 
cio (1450) which would require to include the higher order terms of our model. 
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2 Brief description of the model 

A scalar meson H with the four-momentum P and the mass M is covariantly represented 

MM by 

\H(E, P, M)> = ^ EEC// [4m a m 2 ] rf 4 p 5(p 2 - m?) 6(e) 



• d 4 q 5{q 2 - m 2 ) 0(e) d 4 K F(K) 5 {4 \p + q + K - P) ®{E) <f) f (l ± ) (2.1) 

■^f,sSP) v hM d+(q,c,s 2 )b+(p,c,s 1 )\0) 
Here m ; are quark masses and / stands for quark flavor. The quark wave function is 

MH) = — ^— (2.2) 

with Z M = (p — qY /2, p M = (e,p) and g M = (e, q). The fitting parameter (/ = 

w, d, s, c) is fixed by fitting the meson mass as described below. The dipole form (2.2) 
was found to be a better choice than the exponential form used earlier ]18], [0|. (See also 
some remarks in Appendix.) Coefficient Cf indicates the flavor content of a particular 
meson (For example, see below formula (3.14), where for / (980) C u = cos(9) /\/2). 

The symbols u,v,d + ,b + correspond to valence quarks while the sea function has a 
general form 



F(K) = - ^ PV{P J q) » ] tp(K). (2.3) 



For simplicity we set 



<p(K) = 1. (2.4) 

The complex looking Dirac function in (2.3) simplifies the model structure easing all for- 
mal manipulations. One could produce a somewhat more complicated model, without 
that Dirac function. Additional parameter(s) in the sea model function would lead to a 



richer and more flexible model [19|. Thus the choice (2.4) correspond to a minimalistic 
model version. 

The scalar meson state (2.1) is normalized so that the matrix element of the vector 
current would be, for example, 

mPfWmPi)) = ^F + (Q 2 )(p f + p.y + ■ ■ ■ (2.5) 

with Q M = (Pf — Pj)^, F+{Q 2 = 0) = 1. The normalization (2.5) insures that the vector 
current, and thus charge, is conserved. This requirement is equivalent to the condition 

(H(E, P, M)\H(E, P, M)) =2E = 
_N(P) 2 r mi m 2 M <f> f (l±) ja /™i™2 M<t> r {V ± ) 



l Z7r J c,ai,s 2 c',s' v s' 2 /,/' J e VII 



e> E 



Here 



[Vf>,c>,s> 2 (^)w/',c',4 (P')uf,c, si (P)vf,c,s 2 (q)}- 



{0\b r (p', c', s'^dfiq', c', s' 2 )d}(q, c, s 2 )bj(p, c, Si)|0)|^. =Tii j =T , (2.6) 



P 



^jml-ml+p\ p. p p. q 

T = 1 + n ; n = ^ ; g " = ^~ m 

After a lengthy but straightforward manipulation, one find 
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N(P) = ^iV(O) 



(H(E, P, M)\H(E, P, M)) = 2E = 3N(0) 2 Y C 2 f [ d 3 p - (Mil) 2 ( pq - mi m 2 ) (2.8) 

1 y J e m 

^From (2.8) N(0) can be calculated numerically. 

The matrix element of the conserved vector current V 1 has to vanish when current 
acts on the scalar meson state, i.e. 



(0\V\H(E,P,M)) =0 (2.9) 

The model states (2.1) are consistent with this very general requirement. Some addi- 
tional details are shown in the Appendix. 

So far the model is closely related to ISGW model p0| . In the nonrelativistic limit 



and in the weak binding approximation it goes exactly in the ISGW form [18], |T9|| . 
However the weak binding approximation means that the quark masses and the quark 



energies are approximately equal POfl . In the present application, the model quark 
fields enter a loop (Fig. 1) which constitutes the lowest QED approximation. The 
QCD corrections are modeled by the functions (2.2) and (2.3). One can try to mimick 
the real QED/QCD world by retaining small (current) quark masses in the model. Then 
the meson mass should be equal to a sum (weighted by the sea influence) of average 
model quark energies. In the model determined by (2.2) and (2.4) this i 



M = E Cj I d*p (MO f P ■ g - ^ (p,, + U , (2.10) 

IE f f J q\\ e/e V 11 ^ i|;i 9=-p+^(pii)^ v ' 



3iV(0) 2 n2 f ^ ,<j>f(l±) N2 P ■ q - m x m 2 _ 
/ 

As discussed in the Appendix this simple form holds in the minimalistic model version 
(2.4). The wave function P can be connected with the usual potential 17,19,20 as shown 
in Appendix. In the nonrelativistic, weak binding limit (WBL) (2.10) goes into 



M= ((e + e)) W -^ L m l + m 2 (2.11) 
Here rhi are constituent quark masses [BO]. This WBL makes sense only if one uses 



constituent masses m,, with the corresponding /3's [20] in all relevant formulae. 

If (2.10) is calculated explicitly in our model it can hold only for particular values 
of model parameters, i.e. {3f s H with a particular set of quark masses. When one aims 
for rrti close to the current quark masses, that requires the consistent /3f t H values. As 
explained in Appendix, by using Eq.(2.10), one determines the theoretical form factor 
at the physical momentum transfer Q 2 . However, one is still dealing with some sort of 
a mock meson description. 

Small quark masses, which are used with the relativistic model, lead to the (3 values 
which are close to those used in the nonrelativistic model. Full comparisment between 
those cases is given in Appendix. 
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3 Electromagnetic widths 



In our model [18, [T!| the amplitude Ai for the transition — > 27 is determined from 
the leading diagrams shown in Fig. 1. 





Fig. la) 



Fig. lb) 



Figure 1: Two-photon decay. Full lines are valence quarks, wavy lines are emitted 
photons and blobs symbolize scalar meson states. 



Although the diagrams in Fig. 1 closely resemble the free quark diagrams, they are 
not the same. One has to sum over the moment p, q and over the spins of the valence 
quark states. These sums are weighted by the corresponding functions in the meson 
state (2.1). A loop corresponds to each quark flavor. For example, for the flavor d the 
amplitude corresponding to the diagram in Fig. la is determined by 
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Mr = (2tt) 3 / 2 (0| : VrfSpih - h)Y*d ■ \(dd)) = 



/ct'^li'^ TTIj^ > > > > 

m. : V ) f 



7 M-Jttf-m? 



: ^^[b d (l 2 , (3, c^Udik, (3) + d~^(l 2 , (3, ci)v d (l 2 , (3)} : 



(2tt)3 I 

l Z7r J c,si,s 2 J e y|| 
The contraction of the creation (annihilation) operators in (3.1) leads to the summation 
over spin indicies. That give 



Ir[l (-2ph) 7 2m d 2m d 1 W+£(p„)t " ^ V 1 r ^ ^ 
The amplitude corresponding to Fig. 1 is 

M = MfVfe, A 2 )e,(A; 1 , Ai) + Mfe^h, \i)e v (k 2 , A 2 ), (3.3) 

= I(p,q)Tr(k 1 ^k 2 y 1 ' 
Routine calculation of traces produces the final result which contains 

Tr{hY v = A x (2fq» - 2 V »f + k?(p - q) v + h\(p - q)„ - g^[h(p - q)}) (3.4) 



It is convinient to carry out further calculation in the meson rest frame (MRF). That 
is defined by 

= (M,0) ; e = e ; q = -p 

M 

ki = (to, k) ; k 2 = (u;, —A;) ; u — \k\ — — (3.5) 
Further simplification is obtained by selecting orthogonal polarization vectors 

P • ei )2 = h ■ €1 = k 2 ■ e 2 = 



e"(k 2: A 2 ) = 6% = (0, e 2 ) ; e^, Ai) = e? = (0, ei) ; fc 2 • ei = h • e 2 = (3.6) 

As shown in the next section the result does not depend on a particular gauge. With 
(3.5), (3.6) one obtain 



M = -21 (p, ^){A 1 [(e 1 e 2 )(kp^ + 2(pe 1 )(pe 2 )] + A 2 [ - (m)(kp) + 2(pe 1 )(pe 2 )]} (3.7) 



— * 

By using pA; = pucosO one obtains 

.. 3iV(0) /■ 2 . m 2 , 1 ^^,2p\uJcos 2 d + esin 2 d) 



m d u(27iy/ 2 J " * e 2 (1 + p 2 /(4/3| /o )) 2V 1 ' 7 e 2 - p 2 cos 2 # 



= (eie 2 ) • ^(/o) 
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(3.8) 



The calculation of the decay width 



32tiM h 

requires the summation over photon polarization states, 



\M H \z (3.9) 



El^| 2 = 2-/ />7 (^) 2 (3.10) 

A1A2 

as well as the summation over quark flavors in (3.8), connected with the meson quark 
structure, which we parameterize a \H) = J2fCf\f,f). 

cosO 

|/o(980)> = — t-(|ww) + \dd)) + sin9\ss) 
a/2 



|ao(980)) = —f={\uu) - \dd)) 
v2 



— SITlO 

|/ (1370)> = — i^{\uu) + \dd)) + cos6\ss) 
V2 



| Xco (3415)> = \cc) (3.11) 
The physical mixing for these states is usually determined WM by 



cos6 = l; sin9 = (3.12) 



or 



cosd = - ; sin9 = — ^— (3.13) 
o o 

Eventually one finds by summing over flavors: 

In*(/o(980)) = -[^=-/ u5 (/ (980)) + sin0/ a5 (/ o (98O))] 
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e 2 



Jrat(a (980)) = — / uS (oo(980)) 



7^(/o(1370)) = _[_-^/ un (/ (l370)) + cos9I s - s (f (1370))} 



4e 2 

Jnt( Xco (3415)) = — 7 c5 (xco(3415)) (3.14) 

where I^j{H) is given by formula (3.8) and e is the electron unit charge , i.e. e 2 / (Air) = 
a = (137, 04)- 1 . 

The 27-decay width can be put in the final form valid for a meson H from (3.11): 



ira 2 Int(H). 



m) = M-i—^r (3-15) 
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4 Gauge invariance explicitly tested 

Through its covariant nature, being in a sense a simplified rendering of the real QCD 
field theory, the model ]TJ, [TjJ automatically produces gauge invariant results. 
Under the gauge transformation 

e iM -> + Ak if , (4.1) 
the amplitude (3.3) should not change. That leads to the equality 

Tr^r (e 2 + AA^Uex + Ah) u + Tr(h - fc 2 )^( ei + Afc 1 ) M (e 2 + A£; 2 )„ = 



= Tr{k x y v e 2il e Xv + XV (&i -> fca)"^^ + iV(A, fc 2 ) (4.2) 

Here first two terms give the amplitude (3.3). The piece N(A, ki, k 2 ) should not con- 
tribute to the physical amplitude (3.7). In the meson rest frame, using (3.5) and (3.6), 
one immediately finds 

N(A, k u k 2 ) = A— (-?*L + J*L)[ U - e{p)} (4.3) 
m d p ■ ki p- k 2 

while the terms proportional to A 2 cancel. AT (A, k\, k 2 ) enters the integration over the 
bound quark momentum p, as shown in (3.2), (3.3). The corresponding change in A4 
is 

AM = I(p, q)\ p , ■ AT (A, h, k 2 ). (4.4) 



Here one has 



d 3 p = / p 2 dp / sinOdO / d$ 

jO JO JO 



f>£i,2 = \p\[sin9cosQ{xeifl) + sin9sin§(yei t2 )] (4.5) 
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Integration over the azimuthal angle $ gives zero result, so one has 



AM = (4.6) 

as required by the gauge invariance. 
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5 Results and discussion 



The application of the model (2.1) starts with the self-consistency condition (SCC) 
(2.10). Quark masses are selected, as close as practical to the Particle Data [|T| val- 
ues. Than parameters /3f t u are varied for various flavors appearing in (3.11) so as to 
reproduce the experimental meson masses [jlj. 

M[/ o (980)] = 0.980 GeV ; M[a (980)] = 0.980 GeV 

M[fo(1370)] = 1.370 GeV ; M[x co (3415)] = 3.415 GeV (5.1) 

Theoretical expression (2.10), (2.11) is a sum of parts corresponding to various 
flavors. For example 

M[a (980)] th = i(M M + M d ). (5.2) 

It turns out that SCC requires light quark masses somewhat larger than Particle 
Data pi] median values, while strange and charm masses could be kept within Particle 
Data limit. The satisfactory selection i 

m u = m d = 0.015 GeV m s = 0.120 GeV m c = 1.5 GeV (5.3) 

The most stringent restriction on the light quark parameters is obtained by fitting the 
mass of the a$ meson. If there is no strangeness mixing, i.e. with 9 = in (3.11), the 
same mass is calculated for the /o(980) meson too. Some interesting values are shown 
in Table 1. 

With ideal mixing (8 = 0) /o(1370) has a pure ss configuration. The corresponding 
mass values are shown in Table 2. 

The mass of the pure cc state Xc can be reproduced by using (3 C = 0.267 GeV as 
shown in Table 3. 
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Table 1: Mock mass values for ao 



PuA GeV ) 


Mo (Gel/ j 


0.260 


0.885 


0.270 


0.919 


0.280 


0.953 


0.288 


0.980 


0.295 


1.004 



Table 2: Mock mass values for ss configuration 



(3 s (GeV) 


M {GeV) 


0.350 


1.272 


0.381 


1.370 


0.400 


1.435 


0.450 


1.599 


0.500 


1.764 



When nonideal mixing (3.11), (3.13) is allowed the masses of / (980) and / (1370) 
can be reproduced by f3 u and (5 S which are different from those shown in Tables 1 and 2. 
However, the values in Table 1 still correspond to the a mass. The masses M[/ o (980)] 
and M[/ o (1370)] are reproduced by: 

(3 u4 = 0.419 GeV ; (3 S = 0.242 GeV (5.4) 

The corresponding T(H — > 27) values are summarized in Table 4. 

All conclusions depend strongly on the quark masses. For example, if one chooses 
m c = 1.4 GeV than the SCC requires (3 C = 0.346 GeV. 

The qq structure, which is the main feature of the model, might be capable of 
explaining two photon decay. By that one does not mean a naive " free-quark" structure 
of an early nonrelativistic model. In the present model the valence quarks are immersed 
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Table 3: Mock mass values for cc configuration 



n ( /~i t t\ 

(3 c {GeV) 


M (GeV) 


0.250 


3.377 


0.267 


3.415 


0.280 


3.445 


0.300 


3.490 


0.330 


3.560 



Table 4: Decay widths 



Meson 


Mixing 


Ftheory{keV) 


F exp (keV) 


a (980) 




0.137 


0.26 ±0.08 


/o(980) 


(3.12) 


0.380 


0.56 ±0.11 


/o(980) 


(3.13) 


0.534 


0.56 ±0.11 


/o(1370) 


(3.12) 


0.348 




/o(1370) 


(3.13) 


0.145 




Xc (3415) 




4.608 


4.0 ±2.8 



in a sea, which, however rudimentary, takes into the account the interference of other 
QCD induced configurations like for example ss pairs, gluons etc. The SCC (2.10) 
transmits that into the numerical results. 

The experimental error in / (980) — ► 2 / y rate is rather large. Although, the large 
theoretical prediction in Table 4, seems to be in better agreement with experiments, the 
smaller one, which corresponds to the ideal mixing, cannot be ruled out. However, the 
/o(1370) decay into pions indicates the presence of the light qq combinations^, [II], [EL 



Our result also agrees with the nonideal mixing as considered by Lanikjjn]]. If the 
corresponding theoretical predictions (Table 4) for the decay of /o(1370) turns out to 
be at least approximately correct, one would have a very strong support for nonideal 



mixing |]13, |14[. 



The experimental data |jj for the ao(980) — > 27 decay width contain large errors. 
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Our theoretical value, Table 4, is close to the lower experimental limit. Various other 
theoretical approaches are summarized in Ref.(9). Our approach has some analogy with 
Deakin et al.fjHJ who used constituent quark masses and concluded that theoretical 
results depend strongly on the numerical values of those masses. The same strong 
dependence on the masses, was found here. The decay width and the mass of Xc are 
very well reproduced within the model. 

Naturally all our conclusions depend on the validity of model as such. Here we have 
tried a simplicistic version of the model, which relies on the functions (2.2) and (2.4) 
That gave SCC (2.10) which represents a strong restriction on the model parameters. 
By selecting other functions instead (2.2) and (2.4) one would end with less restrictive 
SCC. 

The model which was employed here indicates the importance of the valence qq 



structure [^0J in the meson state. There is some hope that such relativistic model, at 
least in some richer version, can play a useful role in the classification of meson states. 
It can be a useful tool in the design of future experiments if it can provide a reasonable 
estimate of the magnitude of expected experimental effects. 
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Appendix 



If one evaluates the expression (2.10) for arbitrary (3/s, one obtains a value Mo, which 
is different from the physical mass M. 

Thus the model meson state (2.1) is an approximation of the real physical state 
(of a meson with mass M and momentum P M ) in the sense that there is a one to 
one correspondence between physical state with velocity = P^/M (with respect to 
the meson rest frame) and model state with the same velocity This can be seen 
explicitly from formulae (2.6) and (2.7), where the frame dependence of the internal 
quark momenta is described only through the velocity components E/M and P/M 
and/or through the Lorentz scalar quantities 

Mq has some similarity with so called "mock mass" f2"0fl . Therefore, a model state 



(2.1) correspond also to the different momentum Pq given by 

PS = «" ■ M = — • M . (Al) 

Consequently, when calculating a physical quantity dependent on a square of the 
physical momentum transfer Q 2 , one obtains the value of that quantity at the momen- 
tum transfer Qq which is shifted by the factor Mq/M 2 to the physical one, 

Ql = §- Q 2 - (A2) 

This shift is especially important in processes described by the hadronic matrix 
element of the form (0|r M |y4), as for example in leptonic meson decays or A — > 77 
transitions etc. If M5 2 7^ M 2 one obtains amplitude T(Mq) instead of T(M 2 ). For 
light mesons (tt,K...) approximation of T(M 2 ) by T(Mq) is poor. Thus, Ref.(20) has 
introduced suitable corrections. 

For heavy mesons such approximation is much better and so it was not even men- 
tioned in our previous work [0, pi]. 
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In the nonrelativistic quark model [20] "mock-mass" was defined simply as a sum of 
the constituent quark masses. A formal, covariant expression for that is the expectation 
value of the valence quark (antiquark) momentum operators 

d 3 k rrif 



(2tt) 3 k° 



F [b + f {k 1 s)bf{k ) s) + d+{k 1 s)d f {k,s)\ . 



One obtain 



H = E *? (A3) 



7^(H(E, P, M)\P$\H{E, P, M))\p#i B = ^-M (A4) 



mo)2 ECf [d 3 P (^) 2 P ' q -j nim w+<n\ i 



P» 3iV(0) 2 v 2 f 3 <p f (h) 2 p-q-m im2 



9|| ' e/e ^" ' ^^'^-P+^bii)^ 

Here M is a Lorentz scalar quantity which satisfies (P^) 2 = Mq. However, this does 
not mean that our state (2.1) is really an eigenstate of the meson four-momentum. The 
operator P M , which contains only the free quark operators, is a mock-meson operator. 
Our model is Lorentz covariant but it is not relativistic in the quantum field theory 
sense. The model mocks a hyperplane projected solution of a Bethe-Salpeter equation. 
It corresponds to a quasi potential approximation. 
For a "real" meson, one should have 

^(H(E,P,M)\P^\H(E,P,M))\ 0=PftH = P" (A5) 

Here, the quotation marks symbolize the pseudo realistic ( mock ) character of a meson 
state. In the rest frame P M = (M, 0, 0, 0) this determines a mock mass Mq, which is 
given by (2.10) ( M — > Mq ). The expression (2.10) is a normalization integral (2.8), 
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multiplied by the factor py + qu, which is, in the rest frame, the sum of quark energies 
e + e. 



One reaches WBL by introducing large constituent masses m« [20| and by going into 
nonrelativistic limit: 

Mo = N 2 WBL £ / d 3 p (p) 2 0/(U) 2 (^i + m 2 ) = rh x + m 2 (A6) 
/ 

Working with a more general expression (A4) one can enforce M = M by appro- 
priate choice of the model parameters f3 and mj. That choice of parameters ensures 
that one works at the kinematically correct point Ql = Q 2 . 

The present model parameter (3 is dependent not only on the quark-flavours, but also 
on the meson mass. Thus for example f3 u has different values in /o(980) and /o(1370) 
mesons. 



In the simplest version of the model |L8, [X9|] considered here, the role of the quark- 



gluon sea described by the momentum is mostly kinematical. In principle the sea 
could enter into (2.1) dynamically also, affecting both, the internal momentum distri- 
bution and the internal spin distribution. These possibilities, sketched in Ref.(18), 
are not explored here. To some extent their effects were taken into account phenomeno- 
logically by fitting the parameter 

The model parameters can be also connected with the usual Coulomb plus linear 



potential [Eg, 21 



V(r) = -— + br + c (A7) 
3 r 

The comparison with the pseudoscalar meson applications |2(J will be facilitated if 
that case is briefly revised first. The relativistic case is described by the formulae of 
Ref.(23) in which the potential (|A7|) must be used. Their formulae are pseudoscalar 
version of our expression ( |A8| ) below. 

The relativistic model fit requires slight readjustment of parameters. One has to 
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use c = in order to reproduce masses. However the (3 values, which were found by 
variational procedure[20] do not depend on c. In Fig. A.l. the /3's for the relativistic 
and nonrelativistic[20] fit are compared. 




In the relativistic approach one finds sizable f3 values (comparable with nonrealtivis- 
tic ones), for small quark masses also. For qq pairs the ratios among calculated masses 
are for example m^R^ud) /mR{ud) = 0.7, rriNR{us) /m^us) = 0.8. 

The scalar meson masses can be connected with (IS/?]) by 



EfCj g{P,m q ) 

poo 2p 2 6 2 f 
f(f3,m q ) = ^p 2 dp-f^2e f (A9) 
jo p 2 + mj 

h(P,m q )= W&V o f 2 V{p) (A10) 
Jo p 2 + mj 

g(P,m 9 ) = r Wdp l P ^ f 2 (All) 
jo p 2 + mj 

This is just the expression (A. 4) in meson rest frame with ( |A7| ) included and g is 



normalization. With the potential (|A7|) (c = 0, a = 0.63) one finds for example for 
a (980) : (3 = 0.284 and M = 977 MeV. It is useful to note that one obtains h = 
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for (5 = 0.288 and T,h/g = 0.011 MeV for (5 = 0.284. Thus the expressions (A.4) and 
(A. 8) are numerically consistent for our range of parameters. 
The model also satisfies a very general condition (2.9). Using 



V(0)= £ : ^(0)7^/(0): (A12) 
/ 



where 



*/(«) = E / (^s ! ^[bf(l',s)u f (l, S ) + dj(l, S )v f (l, S )} (A13) 
the condition (2.12) can be explicitly written as: 



{0 „ 0} ^ mql 'M^ (A14) 

The conserved vector current requires equal quark masses, i.e m x = m 2 . In the rest 

— * 

frame E = M, P = 0, p = —q, e = e. For the time component of (A6) one obtains 

miq ° ~ m2P ° = ^ = 0. (A15) 

171x1712 m l 

The spatial components /i = 1, 2, 3 vanish after the integration over d 3 p in I(p, q) (3.7). 
A more general scalar meson state than (2.1) can be constructed by replacing 

P 

UfVf -> uf(a + b—)v f (A16) 

and by adjusting the normalization accordingly. It can be shown that the term pro- 
portional to b neither influences the normalization, nor contributes to the decay width 
(3.15) within conventions, such as for example (2.4), which were used here. Therefore 
the simpler form (2.1) was used. 
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